arXiv:1501.02186vl [cond-mat.soft] 9 Jan 2015 


Asymmetric crystallization during cooling and heating in model glass-forming systems 

Minglei Wang,^’^ Kai Zhang/’^ Zhusong Li/ Yanhui Liu/’^ Jan 
Schroers/’^ Mark D. Shattuck/’^ and Corey S. O’Hern^’^’^’^ 

^Department of Meehanieal Engineering and Materials Seienee, 

Yale University, New Haven, Conneetieut, 06520, USA 
^Center for Researeh on Interfaee Struetures and Phenomena, 

Yale University, New Haven, Conneetieut, 06520, USA 
^Department of Physies and Benjamin Levieh Institute, 

The City College of the City University of New York, New York, New York, 10031, USA 
"^Department of Physies, Yale University, New Haven, Conneetieut, 06520, USA 
^Department of Applied Physies, Yale University, New Haven, Conneetieut, 06520, USA 

(Dated: January 12, 2015) 

We perform molecular dynamics (MD) simulations of the crystallization process in binary 
Lennard-Jones systems during heating and cooling to investigate atomic-scale crystallization ki¬ 
netics in glass-forming materials. For the cooling protocol, we prepared equilibrated liquids above 
the liquidus temperature Ti and cooled each sample to zero temperature at rate Rc. For the heating 
protocol, we first cooled equilibrated liquids to zero temperature at rate Rp and then heated the 
samples to temperature T > at rate Rh- We measured the critical heating and cooling rates R^ 
and R*, below which the systems begin to form a substantial fraction of crystalline clusters during 
the heating and cooling protocols. We show that Rl^ > R*, and that the asymmetry ratio R^/Rc 
includes an intrinsic contribution that increases with the glass-forming ability (GFA) of the system 
and a preparation-rate dependent contribution that increases strongly as Rp R* from above. We 
also show that the predictions from classical nucleation theory (CNT) can qualitatively describe the 
dependence of the asymmetry ratio on the GFA and preparation rate Rp from the MD simulations 
and results for the asymmetry ratio measured in Zr- and Au-based bulk metallic glasses (BMG). 

This work emphasizes the need for and benefits of an improved understanding of crystallization 
processes in BMGs and other glass-forming systems. 

PACS numbers: 64.70.pe,64.70.Q-,61.43.Fs,61.66.Dk 


I. INTRODUCTION 

Crystallization, during which a material transforms 
from a dense, amorphous liquid to a crystalline solid, oc¬ 
curs via the nucleation and subsequent growth of small 
crystalline domains [l|. Crystallization in metals has 
been intensely studied over the past several decades with 
the goal of developing the ability to tune the microstruc¬ 
ture to optimize the mechanical properties of metal al¬ 
loys M- However, in-situ observation of crystallization 
in metallic melts is limited due to the rapid crystalliza¬ 
tion kinetics of metals IfQ. 

In contrast, bulk metallic glasses (BMGs), which are 
amorphous metal alloys, can be supercooled to temper¬ 
atures below the solidus temperature Tg and persist in a 
dense, amorphous liquid state over more than 12 orders of 
magnitude in time scales or viscosity [H . Deep supercool¬ 
ing of BMGs provides the ability to study crystallization 
on time scales that are accessible to experiments [Ml. 

These prior experimental studies have uncovered fun¬ 
damental questions concerning crystallization kinetics in 
BMGs. For example, when a BMG in the glass state 
is heated to a temperature Tf < Tg in the supercooled 
liquid region, crystallization is much faster than crystal¬ 
lization that occurs when the metastable melt is cooled 
to the same temperature Tf [H, [ij • Asymmetries in the 
crystallization time scales upon heating versus cooling of 
up to two orders of magnitude have been reported in ex¬ 


periments [m, [T^. The asymmetry impacts industrial 
applications of BMGs because rapid crystallization upon 
heating limits the ther mop lastic forming processing time 
window for BMGs 

Recent studies have suggested that the asymmetry in 
the crystallization time scales originates from the tem¬ 
perature dependence of the the nucleation and growth 
rates , i.e. that the nucleation rate is maximal at a 
temperature below that at which the growth rate is max¬ 
imal. According to this argument, crystallization upon 
heating is faster because of the growth of the nascent 
crystal nuclei that formed during the thermal quench to 
the glass. In contrast, crystallization is slower upon cool¬ 
ing since crystal nuclei are not able to form at high tem¬ 
peratures in the melt. However, there has been no direct 
visualization of the crystallization process in BMGs, and 
it is not yet understood why the asymmetry varies from 
one BMG to another and how sensitively the asym¬ 
metry depends on the cooling rate Rp used to prepare 
the glass. An improved, predictive understanding of the 
crystallization process in BMGs will aid the design of new 
BMG-forming alloys with small crystallization asymme¬ 
try ratios and large thermoplastic processing time win¬ 
dows. 

We employ molecular dynamics (MD) simulations of 
bidisperse spheres interacting via Lennard-Jones poten¬ 
tials [2214^ to visualize directly the crystallization pro¬ 
cess upon heating and cooling in model metallic glass- 
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FIG. 1: (Color online) Snapshots of the nucleation and growth of crystal clusters at several temperatures T as a monodisperse 
Lennard-Jones system is heated from zero temperature to T/ = 2.0 at a rate Rh < Rh (top row) and cooled from initial 
temperature Ti — 2.0 to zero temperature at a rate Rc < R* (bottom row). Distinct, disconnected crystal nuclei are shaded 
different colors. The far right panel indicates the number of clusters Nc normalized by la\ of as a function of temperature 
during a typical heating (top) and cooling (bottom) trajectory. The maximum number of clusters is indicated by the 

horizontal dashed line. 


forming systems. We perform thermal quenches of the 
system from a high temperature Ti in the equilibrated 
liquid regime to a glass at T/ = 0 and vary the cool¬ 
ing rate Rc by several orders of magnitude. For cooling 
rates below the critical cooling rate Rc < the system 
begins to crystallize, whereas for Rc > i?*, the system 
remains amorphous. We also performed MD simulations 
in which we heat the zero-temperature glassy states (pre¬ 
pared at cooling rate Rp > R"^) through the supercooled 
liquid regime to Tf = Ti over a range of heating rates Rh • 
For heating rates Rh < Rh, the system begins to crys¬ 
tallize, whereas for Rh > Rh, it remains amorphous. We 
also find that the critical heating rate has an intrinsic 
contribution i?^(oc) and an i?p-dependent contribution 
R'^{Rp) — i?^(oc) that increases with decreasing Rp. We 
measured the asymmetry ratio R^^/R* as a function of 
the glass-forming ability (GFA) and Rp for several bi¬ 
nary Lennard-Jones mixtures and find that RJ^/R* > 1 
and the ratio grows with increasing GFA and decreas¬ 
ing Rp. We show that these results are consistent with 
predictions from classical nucleation theory (CNT) that 
the maximal growth rate occurs at a higher temperature 
than the maximal nucleation rate and that the separa¬ 
tion between the nucleation and growth peaks increases 
with the GFA. Further, CNT is able to qualitatively re¬ 
capitulate the dependence of the asymmetry ratio on the 
GFA as measured through R* for both our MD simula¬ 
tions and recent experiments on BMGs as well as on Rp 
for the MD simulations [25[. 

The remainder of the manuscript is organized into 
three sections: Sec. ttH Methods, Sec. HIl Results, and 
Sec.ITVl Conclusion. In Sec. nn we describe the MD sim¬ 
ulations of binary Lennard-Jones mixtures, the compu¬ 
tational methods to detect and structurally characterize 
crystal nuclei, and measurements of the critical cooling 
and heating rates, R* and In Sec. lIIIi we show results 


from MD simulations for the time-temperature transfor¬ 
mation diagram [2^ and the asymmetry ratio R^/R* as 
a function of the glass-forming ability as measured by the 
critical cooling rate R* and the cooling rate used to pre¬ 
pare the zero-temperature glasses Rp. We also compare 
our results for the asymmetry ratio to experimental mea¬ 
surements of the ratio for two BMGs and to predictions 
of the ratio from classical nucleation theory. In Sec. Ml 
we briefly summarize our results and put forward our 
conclusions. 


II. METHODS 

We performed MD simulations of binary Lennard- 
Jones (LJ) mixtures of = Na + Nb spheres with 
mass m at constant volume V = in a. cubic simulation 
box with side length L and periodic boundary conditions. 
We studied mixtures with Na = Nb and diameter ratio 
= (TbIcta < 1- We employed the LJ pairwise interac¬ 
tion potential between spheres i and j: 

uhij) = ieliaij/nA^ - {(Tij/njf], ( 1 ) 

where Vij is their center-to-center separation, e is the 
depth of the minimum in the potential energy u{rij), 
(^ij = fias been truncated and 

shifted so that the potential energy and force vanish for 
separations Vij > 3.5(7^^ [l^. We varied the system 
volume V to fix the packing fraction 0 = 'kcf\{Na + 
a^NB)/QV = 0.5236 Q at each diameter ratio a. For 
most simulations, we considered N = 1372 spheres, but 
we also studied N = 4000 and 8788 to assess finite-size 
effects. Below, energy, length, time, a nd te mperature 
scales are expressed in units of e, aA and e/kB, 

respectively, where the Boltzmann constant kB has been 
set to be unity. 
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A. Cooling and Heating Protocols 

For each particle diameter ratio, which yield differ¬ 
ent glass-forming abilities, we performed MD simula¬ 
tions to cool metastable liquids to zero temperature and 
heat zero-temperature glasses into the metastable liquid 
regime to measure i?* and at which the systems be¬ 
gin to crystallize. To measure i?*, we first equilibrate 
the system at high temperature Ti = 2.0 using a Gaus¬ 
sian constraint thermostat . We then cool the system 
by decreasing the temperature linearly at rate Rc from 
Ti to Tf = 0: 


T(t) = Ti- R,t. (2) 

To measure the critical heating rate Rl^{Rp) at finite 
rate Rp, we first prepare the systems in a glass state 
by cooling them from the high temperature liquid state 
to zero-temperature at rate Rp > i?*. To measure the 
intrinsic critical heating rate i?^(oo), we quench the sys¬ 
tems infinitely fast to zero temperature using conjugate 
gradient energy minimization. For both cases, we heat 
the zero-temperature glasses using a linear ramp 

T{t) = Rht (3) 

until Tf = 2.0. For both heating and cooling protocols, 
we carried out Wot = 1000 independent trajectories and 
averaged the results. 


B. Identification of Crystal Nuclei 


To detect the onset of crystallization in our simula¬ 
tions [m, we differentiate ‘crystal-like’ versus ‘liquid-like’ 
particles based on the value of the area-weighted bond 
orientational order parameter for each particle 11,133. 
We define the complex-valued bond orientational order 
parameter for particle i: 


Qlm ('^) 




2^j=i 


^ij 


( 4 ) 


where Yim(0{fij)^ 0(w)) is the spherical harmonic of de¬ 
gree I and order m, 0{rij) and (l){rij) are the polar and 
azimuthal angles for the vector W , j = 1,..., W gives 
the index of the Voronoi neighbors of particle i, and Aij 
is the area of the face of the Voronoi polyhedron com¬ 
mon to particles i and j. The correlation coefficient [29| 
between the bond orientational order parameters qim{^) 
and qim{j)i where particle j is a Voronoi neighbor of i. 


^ _ X]m=-6 

(EL-6 I P) (EL-6 I P) 

( 5 ) 

is sensitive to face-centered-cubic (FCC) order. When 
Sij > 0.7, i and j are considered ‘connected’. If parti¬ 
cle i has more than 10 connected Voronoi neighbors, it 





FIG. 2: Shifted and normalized probability for crystalliza¬ 
tion {P{Rh,c) — P^c)/{Ph,c ~ P^c) versus the scaled heating 
or cooling rate Circles (squares) in¬ 

dicate data for cooling (heating) for diameter ratios a = 1.0 
(filled symbols) and 0.97 (open symbols). The insets show 
the fraction of crystal-like particles Ncr/N as a function of 
temperature T during cooling (lower left) and heating (up¬ 
per right) for 12 configurations with a = 1.0. The four solid, 
dashed, and dot-dashed curves in each inset correspond to 
cooling and heating trajectories with rates slower than c? 
near R^^c^ and faster than R^^c^ respectively. Trajectories for 
which Ncr/N exceeds 0.5 (above the horizontal dashed line) 
are considered to have crystallized during the heating or cool¬ 
ing protocol. 


is defined as ‘crystal-like’. The ratio Ncr/N gives the 
fraction of crystal-like particles in a given configuration. 
In addition, we also define a crystal cluster as the set of 
crystal-like particles that possess mutual Voronoi neigh¬ 
bors. Distinct crystal clusters that nucleate and grow 
upon heating and cooling are shown in Fig. [H 


C. Probability for Crystallization 


For each diameter ratio and rate, we measure the prob¬ 
ability for crystallization P{Rh^c) = ^x/Wot, where 
Nx is the number of trajectories that crystallized with 
Ncr/N >0.5 during the heating or cooling protocol and 
Wot is the total number of trajectories (c/. insets to 
Fig. [2]). We find that the data for P(W,c) collapses onto 
a sigmoidal scaling function as shown in Fig. [2j 


iPiRn,c) - PZ) _ 1 

pO _ poo 9 

^h,c ^h,c ^ 


1 - tanh I logio ( w 


^1 l^h,c 


R 


'h,c ^ 


(6) 

where P^^ is the probability for crystallization in the 
limit of infinitely fast rates Rh,c co, P^ ^ is the proba- 























4 


.X 10 



FIG. 3: Maximum value of the uumber of crystal clus¬ 

ters Nc{T) uormalized by L^lcr\ (averaged over 1000 tra¬ 
jectories) duriug the cooliug (squares) aud heatiug (circles) 
protocols at rates Rc ~ O.bR* aud Rh ~ O.bR^ for LJ mix¬ 
tures with diameter ratios a = 1.0, 0.97, aud 0.95. For all 
systems, the maximum uumber of crystal clusters is larger for 
the heatiug protocol compared to that for the cooliug proto¬ 
col aud decreases with iucreasiug glass-formiug ability 

(decreasiug a). 


bility for crystallization in the Rh,c 0 limit, is the 

rate at which P{Rh,c) = {Ph,c + '^h,c is the 

stretching factor. We find that k^c ~ 0-25 and Kh ~ 0.2 
for a = 1.0, and these factors increase by only a few per¬ 
cent over the range in a that we consider. We define the 
critical heating and cooling rates and R* by the rates 
at which P{Rh^c) = 0.5, i.e. 


R 


'h,c 


<cio 


K.h,c tanh 


r<c+^^] 

L n,c-^h,c J. 


( 7 ) 


As shown in the insets to Fig. [21 for Rh^c ^ c i^iost of 
the configurations crystallize during heating or cooling. 
In contrast, for Rh^c ^ none of the configurations 

crystallize. 


III. RESULTS 

An advantage of MD simulations is that they can pro¬ 
vide atomic-level structural details of the crystallization 
dynamics that are often difficult to obtain in experi¬ 
ments. In Fig. in we visualize the nucleation and growth 
of clusters of crystal-like particles during the heating and 
cooling simulations. In both cases, the number of clus¬ 
ters reaches a maximum near T ^0.5. In Fig.|3l we show 
the maximum number of clusters (normalized by 

L^la\) that form during the heating and cooling proto¬ 
cols. We find that more crystal clusters form during the 
heating protocol compared to the cooling protocol for all 


particle diameter ratios studied, which is supported by 
the measured time-temperature-transformation (TTT) 
diagram. In addition, we will show below that the asym¬ 
metry ratio R'f^/R* > 1, and that the ratio grows with 
increasing GFA (increasing diameter ratio) and decreas¬ 
ing Rp. We find that CNT can qualitatively describe the 
dependence of the asymmetry ratio on the GFA, as mea¬ 
sured by the critical cooling rate R*, for both our MD 
simulations and recent experiments on BMGs, as well as 
on the preparation cooling rate Rp for the MD simula¬ 
tions. 


A. Intrinsic Asymmetry Ratio 

The critical heating and cooling rates can be obtained 
by fitting the probability for crystallization P{Rh^c) as 
a function of Rh or Rc to the sigmoidal form in Eq. |6l 
We first investigate the minimum value for the asym¬ 
metry ratio i7^(oo)/i7*, which is obtained by taking the 
Rp ^ oo limit. (The asymmetry ratio Rl^{Rp)/R* for fi¬ 
nite preparation rates Rp will be considered in Sec. lIII CT ) 
In Fig.m we plot i7^(oc)/i7* versus R* (for diameter ra¬ 
tios a = 1.0, 0.97, 0.96, 0.95, and 0.93). We find that 
i?*(oo) > i7* for all systems studied, which is consis¬ 
tent with classical nucleation theory (CNT). As shown 
in Fig. [TJ more crystal nuclei form during the heating 
protocol than during the cooling protocol. In addition, 
CNT predicts that the growth rates for crystal nuclei are 
larger during heating compared to cooling. In Sec. IIIIBi 
we will show that both factors contribute to an increased 
probability for crystallization during heating. 

In Fig. m we also show that the asymmetry ratio 
i7^((X))/i7* increases as the critical cooling rate R* de¬ 
creases, or equivalently as the glass-forming ability in¬ 
creases. In the MD simulations, we were able to show 
a correlation between the asymmetry ratio and the crit¬ 
ical cooling rate over roughly an order of magnitude in 
R*. In Sec. iniBl we introduce a model that describes 
qualitatively this dependence of the asymmetry ratio on 
R*. 


B. Classical Nucleation Theory Prediction for the 
Asymmetry Ratio 


In classical nucleation theory (CNT), the formation of 
crystals is a two-step process: 1) fluctuations in the size 
of crystal nuclei that allow them to reach the critical 
radius r* and 2) growth of post-critical nuclei with r > 
r*. To form a critical nucleus, the system must overcome 
a nucleation free energy barrier: 


AG* 


IGtt I]3 
3 AG'2’ 


(8) 


where AG is the bulk Gibbs free energy difference per 
volume (in units of e/cr^) and S is the surface tension 
between the solid and liquid phases (in units of e/a^). 
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FIG. 4: Intrinsic asymmetry ratio versus the crit¬ 

ical cooling rate R* (for diameter ratios a = 1.0, 0.97, 0.96, 
0.95, and 0.93) normalized by = IK/s on a logarithmic 
scale. The inset shows the intrinsic asymmetry ratio versus 
\ogiQRc/Ro on an expanded scale. The filled circles indi¬ 
cate data from the MD simulations and filled squares indicate 
data from experiments on Zr- and Au-based BMGs Em. 
The prediction (Eq.[T2|) from classical nucleation theory (solid 
line) with A' = (87rAT)o)/3a^ = 0.5 (in units of 
E = 0.26, and Qeff = 2.6 interpolates between the MD sim¬ 
ulation data at high R* and experimental data from BMGs 
at low R*. 


FIG. 5: The nucleation I/ADq (solid lines; left axis) and 
growth UalDo (dashed lines; right axis) rates as a function 
of temperature T for increasing values of the glass-forming 
ability (GFA) c = 1.2, 1.1, 1.0, 0.9, 0.8, 0.7, 0.6, and 0.5 (from 
top to bottom) that span the range of diameter ratios from 
a = 1.0 to 0.93. The filled circles indicate the maximum rates 
(/* and U*) for each GFA. As the GFA increases, /* and U* 
decrease and the difference Tu — Tj between the temperatures 
at which the maxima in U{T) and I{T) occur increases (inset). 


We assume that AG = c(T^ — T) where Tm is melt¬ 
ing temperature, — T is the degree of undercooling, 
and c ^ LyjTm is a dimensionless parameter that char¬ 
acterizes the thermodynamic drive to crystallize and will 
be used to tune the GFA of the system (where Ly is the 
latent heat of fusion). Within CNT, the rate of formation 
of critical nuclei (nucleation rate) is given by: 


/ = ADo exp ( — 


Qeff 


exp 


AG^ 


(9) 


where A is an 0(1) constant with units Dq is the 
atomic diffusivity with units a a y^efm^ and Qef / is an ef¬ 
fective activation energy for the diffusivity (with units of 
e). After the nucleation free energy barrier AG* has been 
overcome and crystal nuclei reach r > r*, the growth rate 
of crystal nuclei is given by 


from a = 1.0 to 0.93.) Both /(T) and U{T) are peaked 
with maxima /* and 17* at temperatures Tj and Tu. In 
Fig. [5l we show that as the GFA increases, /* and 7/*, 
as well as T/ and Tu decrease. However, T/ decreases 
faster than T^/, so that the separation between the peaks, 
Tu — Tu increases with GFA. 

To determine the critical heating and cooling rates, 
and i7*, we must calculate the fraction of the samples 
Nx that crystallize and the probability for crystallizing 
P{Rh,c) = where Wot is the total number of 

samples, upon heating and cooling. Within classical nu¬ 
cleation theory, the probability to crystallize upon cool¬ 
ing from Ti to Tf is given by Q : 


U 


exp 





4:71 

Uf 


r / AGvy 

r H T )\ 

, (10) 


J U{T")dT" 


n 3 


dT'. (11) 


where a the characteristic interatomic spacing. 

In Fig. [5l we plot the nucleation IjADo and growth 
rates Ua/Do with Qeff = 2.6 and W ~ 1-40 from MD 
simulations of binary LJ systems [32|, E = 0.26, which is 
typical for BMGs [l5| , while varying the GFA parameter 
from c = 1.2 to 0.5 (corresponding to diameter ratios 


We assume that Ti is above the liquidus temperature 
Tu and Tf is below the glass transition temperature T^, 
where the time required to form crystal nuclei diverges. 
We can rearrange Eq. [11] to solve for the critical cooling 


Ua/Do 
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rate at which P{R*) = 0.5: 


{Kf 



dT' 

AG* 

T' 



1 — exp 


AGV\ 
T" ) 


dT" 


( 12 ) 


where A' = {^T:ADQ)/{Za^) and we assumed that A, 
and a are independent of temperature. A similar expres¬ 
sion for the intrinsic critical heating rate R^{oc) can be 
obtained by reversing the bounds of integration in Eq. [121 

In Fig. m we plot the intrinsic asymmetry ratio 
predicted from Eq.[T2]versus the critical cool¬ 
ing rate R* after choosing the best value A' = 0.5 that 
interpolates between the MD simulation data at high R* 
and experimental data from BMGs at low R*. We find 
that CNT qualitatively captures the increase in the asym¬ 
metry ratio with increasing GFA over a wide range of 
critical cooling rates from IK/s (experiments on BMGs) 
to 10^^K/s (MD simulations of binary LJ systems). A 
comparison of Figs. |4] and [5] reveals that the increase in 
the intrinsic asymmetry ratio is caused by the separation 
of the peaks in the growth and nucleation rates U (T) and 
/(T) that occurs as the GFA increases. Thus, we predict 
an enhanced value for Tjj — Tj in experiments on BMGs 
since the critical cooling rate in experiments is orders of 
magnitude smaller than in the MD simulations. 

The fact that i?^(oc) > R* is also reflected in 

the asymmetry of the “nose” of the time-temperature- 
transformation (TTT) diagram. In Fig |6l we show the 
probability P that the system has crystallized at a given 
temperature T after a waiting time t for monodisperse LJ 
systems. We find that Tmin ^ 0.5-0.6 is the temperature 
at which the waiting time for crystallization is minimized 
and that the time to crystallize is in general longer for 
T < Tmin compared to T > Tmin- Because crystallization 
on average occurs at a higher temperature during heating 
and a lower temperature during cooling, the asymmetry 
in the TTT diagram indicates that slower rates are re¬ 
quired to crystallize during cooling than during heating, 
i.e. R* < Rl- 


C. Asymmetry Ratio for Finite Rp 

In Sec. IIIIBi we assumed that the initial samples {i.e. 
the zero-temperature glasses) for the heating protocol 
were prepared in the Rp ^ oo limit and, thus were purely 
amorphous. How does the asymmetry ratio Rl^{Rp)/R* 
depend on Rp when the preparation cooling rate Rp is 
finite and partial crystalline order can occur in the sam¬ 
ples? In this section, we show results for the asymme¬ 
try ratio Rl^{Rp)/R* for monodisperse systems using a 
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FIG. 6: The time-temperature-transformation (TTT) dia¬ 
gram during cooling is visualized by plotting the probability 
to crystallize P (increasing from light to dark) from 96 sam¬ 
ples as a function of temperature T and waiting time t for 
LJ systems with diameter ratio o = 1.0. A sample is consid¬ 
ered crystalline if the number of crystal-like particles satisfies 
Ncr/N >0.5. The initial states are dense liquids equilibrated 
at T = 2.0. Each initial state is cooled (at rate Rc ^ R*) to 
temperature T < Ti, where Ti ^ 1.4 is the liquidus tempera¬ 
ture, and then run at fixed T for a time t. 


protocol where the samples are quenched from equili¬ 
brated liquid states to zero temperature at a finite rate 
Rp and then heated to temperature Tj at rate R^. (See 
Sec. Ill AL Note that when Rp/R* ~ 1, some of the sam¬ 
ples crystallize during the cooling preparation, yet these 
samples are still included in the calculation of the proba¬ 
bility P{R'^{Rp)) to crystallize. In Fig. (3 we show the re¬ 
sults for the asymmetry ratio R'^{Rp)/R* from MD sim¬ 
ulations. We find that i?^(i?p)/i?* grows rapidly as Rp 
approaches R* from above and reaches a plateau value 
of ^ 1.2 in the limit RpjR^ ^ 1. 

The critical heating rate Rl^{Rp) at finite Rp can also 
be calculated from GNT using an expression similar to 
Eq. [12] with an additional term that accounts for cool¬ 
ing the equilibrated liquid samples to zero temperature 
at a finite rate. In Fig. [T] we show that the asymmetry 
ratio i?^(i?p)/i?* predicted using GNT agrees qualita¬ 
tively with that from the MD simulations. The number 
of crystal nuclei that form during the quench increases 
with decreasing Rp, which causes R'^{Rp)/R* to diverge 
as Rp ^ T*. The predicted intrinsic contribution to the 
asymmetry ratio for Rp ^ R* is small, and R'^{Rp)/R* is 
dominated by the preparation protocol. In contrast, the 
asymmetry ratio i?^(i?p)/i?* ^ 1.2 is dominated by the 
intrinsic contribution in the Rp ^ R* limit. As shown in 
Fig.[4l the size of the intrinsic contribution to the asym¬ 
metry ratio can be tuned by varying the GFA, which con- 
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on BMGs as well as on for the MD simulations. Thus, 
our simulations have addressed how the thermal process¬ 
ing history affects crystallization, which strongly influ¬ 
ences the thermoplastic formability of metallic glasses. 
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FIG. 7: Asymmetry ratio Rl^(Rp)/R* plotted versus the 
preparation cooling rate Rp normalized by the critical cooling 
rate R* from MD simulations with o = 1.0 (filled circles) and 
the prediction from CNT (solid line) with the same parame¬ 
ters used for the fit in Fig. HJand the GFA parameter set to 
c = 1.2. The vertical dashed line indicates Rp = R*. The 
horizontal dashed lines Rl^(Rp)/R* = 1.18 and 1 indicate the 
plateau value in the Rp ^ R* limit and R^ = R*, respec¬ 
tively. The gap between the horizontal dashed and dotted 
lines give the magnitude of the intrinsic asymmetry ratio for 
this particular GFA (c/. Fig. H]). 


trols the separation between the peaks in the nucleation 
/(T) and growth U{T) rates. 


IV. CONCLUSION 


We performed MD simulations of binary Lennard- 
Jones systems to model the crystallization process during 
heating and cooling protocols in metallic glasses. We fo¬ 
cused on measurements of the ratio of the critical heating 
and cooling R* rates, below which crystallization oc¬ 
curs during the heating and cooling trajectories. We find: 
1) > R* for all systems studied, 2) the asymmetry 

ratio R^/R* grows with increasing glass-forming ability 
(GFA), and 3) the critical heating rate R^(Rp) has an in¬ 
trinsic contribution R^(oo) and protocol-dependent con¬ 
tribution R^(Rp) — R^(oc) that increases with decreasing 
cooling rates Rp used to prepare the initial samples at 
zero temperature. We show that these results are consis¬ 
tent with the prediction from classical nucleation theory 
that the maximal growth rate occurs at a higher temper¬ 
ature than the maximal nucleation rate and that the sep¬ 
aration between the peaks in nucleation I{T) and growth 
U{T) rates increases with the GFA. Predictions from 
GNT are able to qualitatively capture the dependence 
of the asymmetry ratio on the GFA as measured through 
R* for both our MD simulations and recent experiments 
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